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SUMMARY 

A rapidly convergent successive approximation process is des- 
cribed that simultaneously determines both composition and tempera- 
ture resulting from a chemical reaction. This method is suitable 
for use with any set of reactants over the c omplete range of mixture 
ratios as long as the products of reaction are ideal gases. An 
approximate treatment of limited amounts of liquids and solids is ■ 
also included. This method is particularly suited to problems having 
a large number of products of reaction and to problems that require 
determination of such properties as specific heat or velocity of 
sound of a dl ssociating mixture . 

The method presented is applicable to a wide variety of problems 
that include (l) combustion at constant pressure or volume; and (2) 
isentropic expansion to an assigned pressure, temperature; or Mach 
number. 


INTRODUCTION 

The theoretical performance of propulsion systems having high 
combustion temperatures can be calculated on the assumption that 
chemical equilibrium exists among the products of reaction. The 
equilibrium composition and the temperature for a system of Np pro- 
ducts of reaction are determined by the simultaneous solution of at 
least Np+1 equations involving dissociation, mass balance, and 
energy or ' entropy balance . This calculation becomes increasingly 
difficult as Np increases. 

The usual method for solving these equations provides a 
successive approximation or trial-and -error process for determining 
the composition at an assumed temperature and pressure. Examples of 
these methods are found in references 1 to 4. When it is desired to 
find the temperature of a system in equilibrium, with a parameter 
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such, as entropy or enthalpy assigned, the composition is usually 
computed at a sequence of temperatures that either converge to the 
correct temperature or are spaced to permit interpolation to obtain 
the correct temperature. 

A rapidly convergent successive approximation process that 
determines composition at an assigned temperature or that s lmult ane- • 
ously determines both composition and temperature for assigned values 
of another parameter, such as enthalpy or entropy, was developed at 
the KACA Lewis laboratory during 1948 and is presented herein. This 
process also permits computation of the partial derivatives required 
to compute such thermodynamic properties as specific heat and 
velocity of sound corresponding to chemical equilibrium. The equa- 
tions are derived that are required for solution of the following 
cases: (l) combustion at constant pressure or volume; and (2) isen- 
tropic expansion to an assigned pressure, temperature, or Mach number 
Examples are given for (1) constant-pressure adiabatic combustion; 

(2) isentropic expansion to an assigned pressure; and (3) isentropio 
expansion to an assigned Mach number. 

This method is particularly suitable for problems having a 
large number of products of reaction and for problems that require 
determination of partial derivatives. Although.it is possible, at 
least in special cases, to devise a procedure that involves less 
numerical computation, the method presented is applicable in a wide 
variety of cases and its numerical application to a given process is 
always simple and essentially the same for all reactions. 


GEHEBAL METHOD 

The thermodynamic state following a specific process, such as 
combustion at constant pressure, can be determined from an appropriate 
combination of the following equations: (a) dissociative equilibrium; 
(b) conservation of mass; (c) conservation of energy; (d) pressure; 
and (e) entropy. Equations (a) and (b) are used to specify chemical 
equilibrium and, when used with any two of the remaining equations, 
define a process. 

The successive approximation procedure presented herein for 
finding the simultaneous solution of a specific combination of the 
aforementioned equations consists of the following steps: 

(a) Estimates of composition and temperature are made and used 
in simple equations to compute the values of error parameters, which 
indicate inconsistency among the estimates of composition and 
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temperature. (These estimates need not "be "based on previous experi- 
ence hut for rapid convergence it is desirable that they he close to 
the final values.) 

(h) A set of linear simultaneous correction equations are given 
that determine a new composition and a new temperature. 

(c) The new composition is used to compute new values of the 
error parameters and step (h) is repeated until the desired accuracy 
is obtained. 

In order to aid in computation, general instructions are given 
that permit construction of the correction equations in matrix form. 
The numerical values of the error parameters are obtained directly 
from the elements of the matrix of coefficients of the unknowns and 
the solution of the matrix equation provides the correction factors 
necessary for determination of the new composition and temperature. 


Equations for Dissociation, Mass, Pressure, and "Volume 

The substances entering a reaction process will be designated 
the reactants and can be represented by the equivalent formula 

z a 0 Y b 0 2 C o ‘ * * 

where the subscripts a Q , b 0 , and c D are proportional to the 
total number of atoms of the elements Z, T, and X, respectively, 
contained in a quantity of the entering substance at the initial con- 
ditions. (A complete list of symbols is included in appendix A.) For 
example, the reactants for a rocket combustion process using 3 moles 
Of ammonia for fuel and 2 moles of nitric acid (HNO3) for 

an oxidant are 


3KH3 + 2HNO3 

and the equivalent formula would be 


H 11 k 5 °6 


where Z, T, and X are the atoms hydrogen, nitrogen, and oxygen, 
respectively, and a Q , b Q , and c 0 are 11, 5, and 6, respectively. 



4 


MCA TN 2113 


The reaction under consideration can "be -written 

A ( z a 0 ^tc?Co) ~ >n l( z ai^hi 2 ci) + n 2( z a2^h2 S C2^ + 


• • • (Z a ^Tb^Zc^) (l) 

where nj_ is the number of moles of the ith molecule or atom. The 
subscripts aj, bi, and c-j_, which can take on only integral values 
or zero, denote the n umb er of Z, Y, and X atoms in the ith mole-, 
cule. For example, if Z, Y, and X again represent hydrogen, 
nitrogen, and oxygen, respectively, the values of a^, bj_, and c^ 
for a water molecule H2O would be 2, 0, and 1, respectively. 
Assumptions are made that the produ ct s of reaction are contained by 
a volume Y numerically equal to the gas constant R times the 
absolute temperature T so that for ideal gases 


Pi = ni 

During the solution of the problem, determination of the number of 
formula weights of the reactants A that are required to balance 
the reaction given by equation (l) is necessary. Products of reac- 
'tion In the gas phase are assumed to be ideal gases that form Ideal 
mixtures and each condensed phase is assumed to have a partial 
p ressure of zero, even when finely divided and suspended In the gas. 
For solids and liquids therefore 


Pi = 0 

As an approximation, the following assumptions are also made: Each 

condensed product is Insoluble in all others, the fugacity of each 
condensed phase Is equal to 1; the total volume occupied by the 
liquids and solids is negligible with respect to the volume occupied 
by the gases; and the liquid and solid particles have the same 
temperature and flow velocity as the gases. 

Dissocation equations . - The equation for the dissociative 
equilibrium among gaseous atoms and molecules can be written as 

a i z + + c-jX Z ai Y- bl X Ci 


( 2 ) 
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and the corresponding equation for the equilibrium constant 
gaseous molecules is 


% = 


P 


i 


Pz ai Py bi Pl Ci 


Ki of 

(3) 


For liquid or solid molecules, assuming the fugacity of each 
condensed phase is equal to 1, 


E i 


_ M *>i °i 
Pz, • Py Pi 


(4) 


where p^, py, and py are the partial pressures of the ' Z, Y, 

and X atoms in equation (l), respectively. The equilibrium con- 
stants can also be expressed in terms of the free-energy 
changes (AF^) i across the dissociation reactions represented by 
equation (2) or 

lo Ss % = (5) 


Beoause the assumed composition may not correspond to that at 
chemical equilibrium, variables k. are conveniently defined so 
that for gaseous molecules 


loge ^ 


logg Pi - ai loge % " *i lo 8e Py - c i lo Se Py 



( 6 ). 


and for liquid or solid molecules 


logg ^ ^ logg P2 t bi log e Py - Ci logg Py- 



(7) 


The value of each must approach 1 as the solution to the 

problem is found. Applying equation (6) or (7) to each molecule 
results in equations, where denotes the number of different- 

types of molecule considered. 

For simplicity of nomenclature and presentation, the equations 
for dissociation are expressed in terms of the atonic gas, which 
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leads to a simple computational procedure when the atomic gases of 
most elements are present in significant quantities. If the atomic 
gas of an element is not to he considered as one of the products of 
reaction, however, a molecule containing the element may he sub- 
stituted in equation (2) in place of the atomic gas. For example, 
carbon monoxide can he used in- place of carbon gas and would then 
be treated as an atom. 

Mass-balance equations . - A mass -balance equation stating the 
conservation of atomic type can be written for each chemical element 
present. Because a composition Is initially assumed, it is con- 
venient to define parameters a, b, c,... as representing the 
summation over the products of reaction of each atomic type per 
equivalent formula . 


I £ a i n i 

(8) 

i £^i^i 

O) 

i 



(10) 


As the solution to the problem is found by successive adjustments of 
the Initial assumptions, the values of a, b, and c approach the 
values of a Q , b Q , and c Q , respectively. The mass-balance equa- 
tions result in N e equations, where N e denotes the number of 
chemical elements . 


Total -pressure equation . - The total pressure P is the sum of 
the partial pressures 


p=-£Pi (n) 

i 

For a process with an assigned pressure, the value of P must 
approach the assigned value P 0 as the solution of the problem is 
found. 


Constant volume . - For processes that occur at constant volume, 
the density of the mixture is constant. The density p is defined 
as 


P 


AMp _ AM r 

~ ~ hF 


( 


( 12 ) 
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where is the molecular weight of the equivalent formula. For 
a reaction process with an assigned density, the value of p must 
approach the assigned value . p D as the solution of the problem is 
found. 


Combustion at Constant Pressure 

For given initial conditions, the temperature and the composition 
following a combustion process are to be found. When chemical energy 
is included in the enthalpy of each substance, the enthalpy of the 
products of reaction following an adiabatic combustion must be equal 
to the enthalpy of the reactants at the initial conditions. An 
arbitrary base may be adopted for assigning absolute values to the 
enthalpy of various substances because only differences are measur- 
able. One such base (reference 4) assigned a chemical energy of zero 
to the oxidized form of the elements at 100° F. For calculations 
with the elements boron, hydrogen, fluorine, oxygen, nitrogen, and 
carbon, assigning a chemical energy of zero to water, oxygen, hydrogen 
fluoride, boron trifluoride, nitrogen, and carbon dioxide at the 
absolute temperature of 0° K was found to be more convenient. By 
using this base, the chemical energies of all other molecular types 
entering a combustion process containing these elements are positive, 
thus avoiding a possible souroe of difficulty that might occur, in 
the recommended method of adjustment when a logarithm of a negative 
number (or zero) might be required. 

Enthalpy of fuel and oxidant. - The enthalpy at initial condi- 
tions of the amount of fuel and oxidant corresponding to the equivalent 

formula Z a Y-u X c is denoted by h Q and is given by the expression 
tt o u o ^o 

ho = nj> (Brp)f + ng (Hp)g (13) 

where n^ and ng are the number of moles of fuel and oxidant, 
respectively, corresponding to the equivalent formula Za 0 Yb 0 Xc 0 and 
(B^) f . and (B^) g are the molar enthalpies of the fuel and the 
oxidant, respectively, at the initial conditions. The molar enthalpy 
Bp is defined by tbe equation 

B? - Pc* ar + Hg 
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where C° is the molar specific heat at constant pressure, and 
is the chemical energy of the substance at an absolute temperature 
of 0° E. 

Enthalpy of products of reaction . - The enthalpy of the products 
of reaction per equivalent formula can be conveniently represented by 
a variable h that is given by the equation 

k = ( 14 ) 

as long as the kinetic energy is negligible. When enthalpy is 
assigned, the value of h approaches h Q as the solution of the 
problem is found by successive adjustments of the estimated quantities. 
If heat were lost, the value of h Q would be accordingly reduced. 

Equations for constant -pres sure combustion . - The equations 
defining the constant-pressure combustion are: 

Type Humber of equations 

Dissociative equilibrium 
Conservation of mass N e 

Constant pressure 1 

Conservation of energy 1 

These equations are to be solved simultaneously. Values for n^, A, 
and T can be estimated and the values of the parameters k^, a, b, 
c, . . . , P, and h can be computed from equations (6) to (11) and 
(14). Corrections are then required for ni. A, and T. 

Correction equations . - The adjustments to n^. A, and T are 
■hwAp by means of a set of Hp+2 correction equations derived from 

equations (6) to (11) and (14) that adjust the estimates by the Newt on - 
Baphson method for solving simultaneous equations (reference 5) . This 
method can be illustrated by a simple example . If Qq_ and Qg are 
functions of q and r, 

Qq_ = fi (<l,r) 


Qg = f 2 (q,r) 
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By taking estimated values, for example q Q and r 0 , each function 
may be expanded in a Taylor’s series about the point (q 0 ,r 0 ) and 
when derivatives of higher order than the first are neglected 


AQ, = ^ Aq + ^ Ar 
dq Sr , 


SQg ^Qg 

A% - -aj A1 + If 4r 

The desired changes AQq_ and AQg can be computed and if the 
partial derivatives can be numerically evaluated, solving for the 
approximate changes in q and r to effect simultaneously the 
desired changes in both Qy and Qg is comparatively simple 
because the equations are linear. 

The equations for log e k^, a, b, c, . . h, and P can 
similarily be expanded in a Taylor’s series; for example, equa- 
tion (6) for log e kj_ for gases expanded in the logarithmic form 
would be 

A log e ki = A logg Pi - a ± A log e ]% - bi A logg Py - A log e p x 



A logg T 


(15) 


The term (AB^/BT)i appears because 

*(#), 


f AHp' 


S logg T " \HP J ± 


where (AE^)^ is determined across the dissociation -reaction equa- 
tion (2). Expanding equation (8) gives 




l i n i 


AA 


(16) 
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Some of these equations are expressed in terms of A logg and 
Ap^, whereas the others are expressed in terms of An^. In order 
to convert to a common variable A log e n^, the following sub- 
stitutions can be made : 

For gases, p.^ = n^ therefore 

A logg Pj^ = A logg nj_ ‘ (17) 

A Taylor’ s expansion of the logarithm of a variable logg q, by 
dropping terms of higher order than the first, yields 

Aq = q A logg q (18) 

The expansion of equations (6) to (ll) and (14) can be written in 
logarithmic variables with the aid of equations (17) and (18), as 
follows : 

For gaseous products, 

A log n^ - a^ A log - bj_ A log n y " c i A lo 8 n x " 

(AH §) 1 

A log T = - log fc, (19) 

ET 1 

log 

( 20 ) 

(21a) 


For liquid or solid products, 

(AHg)i 

-aj^ A log n^ - 'b 1 A log ny - c^ A log n x ^ — A log T = 

For all products of reaction, 

a Q 

XU a. n. A log n, - Aa A log A = Aa log — 

■j_ 1 1 1 a 

X2^i n i A log n^ - Ab A log A = Ab log 


(21b) 
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ScjEi^ A log n-j_ - Ac A log A = Ac log 

(21c) 

- P 0 

2 _ A log n^ = P log — 

(22) 

A log n^ - Ah A log A + A log T 

h„ 

= Ah log -2. 

(23) 

h 

where the desired values of the parameters (log = o, a 

ii 

P> 

o 

• 

are substituted in the form 

A log e k i = ' lo S e k ± 

(24) 

Aa = a log e 

(25) 


and logarithms to the hase 10 are used for convenience. 

Matrix . - The values of A log A log A, and A log T are 

then computed "by means of the matrix shown in figure 1(a). These 
corrections are applied to the initial set of estimates of n^, A, 
and T and the process is repeated until all the given conditions 
are simultaneously satisfied. .This matrix, however, merely provides 
a convenient soheme for solving the simultaneous equations hut any 
of the well-known methods for solving simultaneous equations may he 
used. 


In order to permit rapid solution of the matrix, using the 
arrangement of rows and columns described in appendix B is 
desirable. With this arrangement, a single step reduces the order 
of the matrix by the number of gaseous molecular types. 


Combustion at Constant Volume 

The procedure given for finding the composition and the tem- 
perature of a combustion process at constant pressure can be applied 
to combustion at constant volume with the following changes: 
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(a) The correction equation for pressure is replaced by a 
correction equation for density obtained from equation (12) 

A log A - A log T = log (26) 

(b) The correction equation for conservation of energy must 
be witten in terms of internal energy E^ and thus becomes 

^(E^)^ A log nj_ - Ae A log A + T ^ ( c °)i n i A lo 8 T 

0 

= Ae log — (27) 

e 

where e Q is the assigned internal energy per equivalent formula 
at initial given conditions and C° is the molar specific heat at 

constant volume. Substitution of these two equations in the matrix 
of figure 1(a) will permit the composition anfl the temperature to 
be found for assigned values of density and internal energy. This 
application of constant-volume combustion, which, for e xamp le-, is 
involved in reciprocating engines and pulse -jet engines, has not been 
made at the Lewis laboratory. 


Isentropic Expansion to Assigned Pressure or Temperature 

Assigned pressure . - The calculation of temperature and equili- 
brium composition of the products of reaction following isentropic 
expansion to a fixed pressure involves the simultaneous solution of 
dissociation, conservatlon-of -mass, pressure, and entropy-balance 
equations . 


For the reaction of equation (1), the dissociation, conservation 
of mass, and pressure equations ((6) to (11)) can again be applied. 
For the conditions following an isentropic expansion, the entropy s 
the products of combustion per equivalent formula after expansion 
must be equal to the entropy s Q of the products of combustion per 
equivalent formula before expansion. 


s 


o 


r* 

‘ - BPi lo Se Pi 


(28) 

combustion 

conditions 
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■where (S^) ^ is the absolute entropy of the product i at stand- 
ard conditions. This formula is applicable for ideal solids and 
liquids, assuming p.^ = 0, as long as their -volume is negligible. 

After the expansion takes place, the entropy per equivalent formula 
is given by the expression 



n i (S,9) i - Rp. logg p y (29) 

— rJ exit conditions 


Whereas equation (28) is, of course, evaluated at combustion-chamber 
temperature an 6 pressure, equation (29) is evaluated for exit tempera- 
ture and pressure. As the solution of the problem is found by succes- 
sive adjustment of estimated quantities, the value of s approaches 

8 o* 

In the adjustment of the values of n^. A, and T, the 
correction equations (19) to (22), whioh have been derived from 
equations (6) to (11), can be applied. In addition, the following 
correction equation for entropy can be written from equation (29) : 


V' 


Sb 1 * a log n^ - As A log A +^(C°) i n i A log T = As log 


o 

8 


(30) 


where - , 

B i = “i " Rp i t 1 + log e p i^ 


The values of A log n^, A log A, and A log T are then 

computed by means of the matrix qf figure 1(b), which is identical 
to figure 1(a) except that equation (30) is substituted in the 
last row in place of equation (23). 

Assigned temperature. - For the computation of data for enthalpy 
entropy diagrams and ^or other practical computations, it is often 
necessary to find the exit pressure and composition as a function of 
exit temperature. The procedure required is the same as that des- 
cribed for isentropic expansion to an assigned pressure except that 
the pressure equation and the temperature column are dropped from 
the matrix of figure 1(b) . 
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Isentropic Expansion to Local Telocity of Sound 

The theoretical velocity of sound that includes the effect of 
dissociation can he computed at any point in a nozzle with a modifi- 
cation of the matrix previously derived to obtain the correction 
quantities . 

Telocity of sound . - The velocity of sound u can he defined 
as 

u2 “© s (31) 

where the subscript s denotes the condition of constant entropy. 
The total differential of pressure dP can he found from equa- 
tion (11) 


dP =^dp 
i 


i 


and the total differential of density dp can he found from 
equation (12 ) . 


dp 



ET d 


(32) 


(33) 


Thus, equation (31) can he written 



and by dividing the numerator and the denominator by dT and 
changing to logarithmic variables 


2 


u 


y a i 08e Pj_ 

i p l d logg T 

; AM r f d l°g e A \ 
KT \d logg T ~ J 


s 


(34) 


This expression will permit evaluation of. u^, provided the values 


of the derivatives 


d loge Pj\ 

,d lo^ T ) 


and 


d logg A 
^d logg T 




are found for 
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conditions of chemical equilibrium and for an isentropic process. 
The conditions of chemical equilibrium and constant entropy are 
introduced by -writing the total differentials of equations (6) 
to (10) and (29). The total differential of these equations 
expressed in logarithmic variables and divided by d log e T 
can be written 

For gaseous products. 


d log pj_ d log pg d log py d log p^ 

d log T ~ a i d log T ” ^ d log T ” Ci d log T 


(AHip)^ d log 
' ET d log T 

For liquid and solid products, 

d log P7 d log py d log p v- 

1 d log T 1 d log T 1 d log T 

For all products of reaction. 


(35) 


(AB^) i d log 
ET d log T 


E< 


i n i 


d 

log 

*1 

- Aa 

d 

log 

A 

= Aa 

d 


a 

d 

log 

T 

d 

log 

T 

d 

log 

T 

d 

log 

“i 

- Ab 

d 

log 

A 

s Ab 

d 

i2SL 

b 

d" 

log 

T 

d 

log 

T 


d 

log 

T 

d 

log 

% 

- Ac 

d 

iog 

A 

= Ac 

d 

log 

c 

aT 

log 

T 

d 

log 

T 

d 

log 

T 


d lo 8 ^ d log A y { 0\ _ d log s 

i a i d log T “ As d log T + V 'S i i As d log T 


(36) 

(37a) 
(37b) 
• (37c) 

(38) 


If d log s is taken as 0, s is a constant; if d log a, 
d log b, d log c, and d log k^ are taken as 0, mss is constant, 

atomic types are conserved, and rate of change in composition 
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corresponds to constant values of log k^. With these assumptions, 

equations (35) to (38) constitute a set of simultaneous equations 
for partial derivatives. The augmented matrix formed from these 
equations iB shewn in figure 2 with the sign reversed. The matrix 
is identical to a portion of the matrix of figure 1(h). When 


d log n^ ^ 

d log T J 3 


and 



are determined by means of the 


matrix shown in figure 2, the velocity of sound can he calculated 
from equation (34) . This equation can he applied to mixtures of 
liquid and solid products in equilibrium as long as their volume 
is negligible compared with the volume of the gas mixture and pro- 
vided the liquid and solid particles move in velocity and tempera- 
ture equilibrium with the gas. 


Specific heat. - The specific heat at constant pressure of a 
mixture in equilibrium may be found, from equation (14) as follows: 
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equal to the enthalpy change from the point of negligible kinetic 
energy to the point in question. From this relation, the velocity 
squared at any point 2 is 



2 (H^)jn^ 


(41) 


■where the subscript 2 indicates that the variables are evaluated 
at point 2 in the nozzle. The Mach number M is 


M = I (42) 

u 

For convenience, a parameter h* is defined as representing the sum 
of heat energy plus kinetic energy of the products of reaction per 
equivalent formula : 


(43) 


s 

As the solution of the problem is found by successive adjustments 
of the estimated quantities, h* approaches ho. 

If equation (43) can be expanded in a manner similar to that 
used to obtain equation (23) and if the differentials of derivatives 
are assumed to be negligible, "the correction equation is 


h* = 


2(H^)^n^ 


+ M 




d log e n i 
d log T 


2A 


d logg A 
d log e T 


- I 1 


5Hh, ' A log n, - Ah* A log A + tX 1(C°)-< ’ A log T 
i 1 i v 


where 


V 


"ifel + 


KT 


d log 


d log T 


log A 
log T 


-1’ 


Ah 


T O 

l0g T[* 


(44) 



18 


MCA TN 2113 


<°SV 


(°$)i 


I^Epj! 


4 log 
log T 


2 /41o^A _i\ 
V4 log T / 


Equation (44), together with equations (19) to (21) and (30), 
constitute the correction equations for the isentropic expansion to 
an assigned Mach number. The coefficients of these equations form 
the matrix shown in figure 3. 


In order to carry, out the numerical co mpu tations, values of 
n^. A, and T are estimated far the assigned conditions; the 


> 

value of 


4 lOg 

d log T J B 


and 


d log a\ 
d log Tj s 


are obtained from the matrix 


of figure 2. The numerical values of the elements of the last column 
and bottom row of the matrix of figure 3 are then computed and by 
using the auxiliary matrix already completed for figure 2, the matrix ' 
of figure 3 may be reduoed and the value of the corrections to n^, 

A, and T found. This process can be repeated until the assigned 
conditions are satisfied. 


Throat area of supersonic nozzle . - The process of isentropic 
expansion to a local Mach n um b er of 1 is particularly interesting 
in the determination of the throat area of a nozzle having greater 
t han critical pressure ratio. By assuming that the flow is isen- 
tropio and that chemical equilibrium is maintained throughout the 
expansion process, the flow velocity v at the throat must be 
equal to the velocity of sound u at the throat. The values n^, 
A, T, and u can he found for a Mach number of 1 by use of the 
procedure given. 


The throat area t can be calculated from the equation 

t = BT 
- m AMyU 


(45) 


where m is the mass flow per second. This equation can be applied 
to mixtures of liquid or solid phases in equilibrium provided that 
the volume occupied by the liquid and the solid phases is negligible 
compared with that of the gas phase and that the particles of liquid 
and solid are in thermal and velocity equilibrium with the gas phase. 
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EXAMELE OF COMBUSTION OF DIBORANE WITH FLUORINE OXIDE 


The calculation, of equilibrium temperature and composition of 
the reaction of 1 mole of diborane (BgHg) with 5 moles of fluorine 
oxide (FgO) is illustrated in this example for processes of 

(a) constant -pressure adiabatic combustion 

(b) isentropio expansion to 1 atmosphere 

(c) isentropic expansion to the local velocity of sound 
An equivalent formula of these reactants is 


Z Y -h W. “ Sr F, n Oc- 
a o B o c o d o 2 6 10 5 


■and a Q = 2, = 6, c Q = 10, and d Q = 5, 


The following gaseous products will be considered as the produo ts 
of reaction: boron trifluoride BF3, boron triaxide BgOj, boron 

fluoride BF, boron hydride BH, boron oxide BO, diatomic boron 
Bg, hydrogen Hg, water vapor HgO, hydroxyl radical OH, hydro- 
gen fluoride HF, oxygen Og, fluorine Fg, atomic hydrogen H, 
atomic boron B, atomic fluorine F, and atomic oxygen 0. Ho 
liquids or solids are included. If the products are numbered in the 
order given, they can be identified in the terminology of equation (l) 
as follows: 


BF 3 =» Bq Hq F 3 0 Q 


and therefore 


a^L = 1, bj_ = 0, c^ = 3, and d^ = 0 


Similarly, 


b 2°3 = ®2 So P o °3 


and 

Eig = 2, b 2 = 0, c 2 = 0, and dg = 3 
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All -values of a^, b^, c^, and d^ for this problem, together 

with the thermodynamic properties used, are listed in table I. 
Although these thermodynamic values have since been revised, they 
are adequate for the purpose of this example . 

The enthalpies of liquid BgB^ at 298.16° K and of liquid 
F 2 0 at 128.3° K, the assumed initial conditions, were computed 
from the heats of formation, the heats of vaporization, and the 
heat changes due to temperature changes together with the arbitrary 
chemical energy assigned to the elements. The value of the heats 
of formation were taken' as -44 kilocalories per mole for B 2 H 5 

(reference 6 ) and 5.5 kilocalories per mole for F 2 0 (reference 7). 

(Heat liberated is considered to be negative.) The enthalpy values 
using the base previously described are 

(S 0 29e.l6> liquid BgH 6 - 570 - 149 m00al0rleB * er mole 


( h ° 128.3) liquid F 2 0 = 67.077 kilocalories per mole 

The enthalpy of the amount of fuel and oxidant at initial conditions 
corresponding to the equivalent formula is, from equation (13), 


h 


o 


570.149 + 5 (67.077) = 905.534 


kilocalories 
equivalent formula 


(46) 


The values of b^, c i , d^, and hg cure constant for all 

parts of this example. 


Combustion Process 

The adiabatic combustion process was assumed to occur at a 
constant pressure of 20.4 atmospheres. 

First estimate . - From previous computations or from simple 
calculations using equilibrium constants, estimating reasonable 
values for the composition and the temperature is usually possible. 
This proce dur e is recommended inasmuch as close estimates reduce 
the number of trials that must be made. In order to show that an 
arbitrary- composition that is not based on probable final values of 
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the composition can he used, however, the first estimates for this 
example for n^ and A have been taken equal to 1 mole end a tem- 
perature of 4000° K. The possibility of divergence is discussed in a 
later section. All estimated quantities will be used with three dec- 
imal places to distinguish them from number's that are always integers. 

Correction equations. - The total pressure from equation (11) 
can be computed because p^ = n^, hence 

P = 21 Pi = 16.000 (47) 


The correction equation for pressure (equation (22)) becomes 

1.000 A log 1.000 A log p^ q t 1.000 A log P^p + . . . + 

3 .23 

1.000 A log p 0 = 16.000 log (48) 


The total amount of boron in the products of reaction can be 
determined from equation (8) 

Aa = S a i n i = 1 ^ + 2 n ^ + 11^ + 1 n^ + l 1 ^ + 2:^ + 

lug = 9.000 (49) 

and therefore the number of estimated boron atoms per equivalent 
formula is 


_ 9.0°° 
a = 1.000 


9.000 


because A has been assumed to be 1.000 for the first estimate. 
The coefficient of each term in equation (8) is equal to the num- 
ber of boron atoms in the molecule; the terms that do not contain 
boron atoms are 0. Each term of equation (49) then becomes the 
coefficient for its respective correction term in equation (21a): 


1.000 A log rLgj, + 2.000 A log q + 1.000 A log n^ + 
3 , 2 3 

1.000 A log n^g + 1.000 A log n-gQ + 2.000 A log ng^ + 


1.000 A log n_ - 9.000 A log A = 9.000 -log 

ij * y • 


000 


(50) 
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Similarly, tlie correction equation based on the conservation of 
hydrogen atoms is 

1.000 A log rtgg + 2.000 A log + 2.000 A log 

1.000 A log n^ + 1.000 A log n^ + 1.000 A log n^ - 8.000 A log A 

= 8.000 log — § — .7* (51) 

8.000 v ' 

the conservation of fluorine atoms is 

3.000 A log ngj, + 1.000 A log n^ + 1.000 A log n^ + 

3 

2.000 A log n^ + 1.000 A log ^ - 8.000 A log A = 8.000 log — 12- 

(52) 

and the conservation of oxygen atoms is 

3.000 A log QggOg + 1*000 A lo 8 “bq + 1 ' 00 ° A l0g n H 2 0 + 

1.000 A log n QH + 2.000 A log n^ + 1.000 A log n Q - 9.000 A log A 

= 9.000 A log _JL_ (53) 

9.000 

For dissociation, the numerical value of log can be computed 

directly from equation (6) vith the data of table I and (by using 
logarithms to the base 10 for convenience and by remembering that 
(-AF§/Er)i log e = log K ± ) is 

log k B F3 - lo S ^BF “ log Pb " 3 lo 6 % " 5 ‘ 695 
3 

As the partial pressures of all the constituents have been esti- 
mated to be 1.000 atmosphere, 

lo S *bF 3 = -5.695 
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The corresponding correction equation from, equation (19) is 

A log - A log Pg - 3 A log p^, - (-62.075) A log T = 5.695 

(54) 

In a similar manner, equation (19) for BgOj is 

A log p - 2A log pl - 3A log p n - (-80.593) A log T 
b 2°3 B 

<= - -log kg o = 5.109 
2 3 

For HgO , equation, (19) is 

A log P H ^ 0 - 2A log Pg - A log P 0 - (-29.209) A log T = -log k^Q 

= - 0.347 

Correction equations for dissociation similar to those given 
for BF^, b 2°3^ and ^2° molecules can he -written for each molec- 
ular constituent considered in the reaction (a total of 12 in this 
example). The sum of the enthalpies of the products of reaction, 
as given in table I and determined by equation (14), is 2734.615 
because n^ = 1.000. The heat-balance equation, as given by 

equation (23), is 

(72.172) A log p__ + (233.435) A log pn + . . . + 

B 2 °3 

(79.493) A log p 0 - (2734.615) A log A + (644.651) A log T 


= (2734.615) log 


905.554 

2734.615 


(55) 
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A matrix (fig. 4) can now be constructed with the 12 equations 
similar to equation (54) and with equations (48), (50) to (53), and 
(55). The recommended arrangement and reduction of the matrix are 
given in appendix B. The solution to the correction equations is 
found to be 


4 108 “BFj 


0.708 

4 108 n 0H = 

0.490 

A log n 

=3 

-1.098 

A log n = 

1.378 

b 2°3 



HF 


A logn^ 

~ 

1.116 

A log n A = 
°2 

0.172 

A log n 

BH 

32 

1.664 

A log n = 
^2 

-2.038 

A log n 

B0 

= 

0.613 

A log n = 
H 

1.299 

A log n^ 

rz 

-2.138 

A log ng <= 

0.929 

A log n^ 

zz 

0.014 

A log Up = 

1.221 

4 108 “HgO 

= 

-0.799 

A log n Q = 

1.457 

A log A 

= 

0.123 

A log T = 

0.154 


These values are to be applied to the initial estimates for 
and T according to the equation 

' (log Vsecond = (log n i ) first + A log n i 
estimate estimate 

For example, the second estimate of n^, would be 

3 


n ] = log 1.000 + 0.708 

■“3/ second 
est im ate 



Hi, A, 
(56) 



second 

estimate 


5.105 
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Another set of correction equations "based on the second estimates 
of 'n^, A, and £ are set up and solved "by means of a second 
matrix. The process is repeated until the values of equations (48) 
and (50) to (55) approach 0. For this example, six approximations 
■we're required to give the following final values of n^. A, and T: 


p bf 3 

= 2.6593 

Pqj = 0.6785 

p b 2 o 3 

= 0.1235 

%F = 7-1456 

P BF 

= 0.1936 

p n = 0.9210 
u 2 

P BH 

= 0.0001 

= 0.0003 

P B0 

= 0.1669 

P H = 1.7694 

p u 

B 2 

= 0 

P B = 0.0577 

% 

= 0.1271 

P F = 1.3043 

p h 2 o 

= 0.0627 

Pq = 5 • 1903 

A 

= 1.6622 

T = 4775.5° E 


Discussion of Convergence . - In order to demonstrate the con- 
vergence of the process with large errors in the first estimate, 
the example of the combustion of diborane and fluorine oxide was 
solved by using 1 mole of each product, a value of 1 for A, and 
a temperature of 4000° E for the first estimate. Because these 
first estimates were made without regard for the probable final 
values, large errors were present in the second approximation 
and six approximations were required to eliminate the error. 

These values of parameters a, b, c, d, P, h, and 5 are shown 
in the following table where & is defined as 


8 = log 


. a o 
log =- 


+ log 


o 

b - 


+|l°g ~ 


!og J-I + 


. p o . *o 
log g- + log g- 
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besoms op appeqhmatiohs 


Param- 

eter 

First 

estimate 

Trial number 

Desired 

value 

i 

2 

3 

4 

5 

6 

a 

8 

36.840 

7.005 

6.286 

6.079 

6.002 

6.000 

6.000 


9 

23.346 

11.605 

2.653 

2.325 

2.008 

2.000 

2.000 

c 

8 

51.540 

24.082 

13.104 

10.541 

10.016 

10.000 

10.000 

d 

9 

29.641 

11.954 

33.660 

5.240 

5.022 

5.000 

5.000 

P 

16 

125.485 

-38.000 

52.434 

21.416 

20.436 

20. 400 

20.400 

li 

2734.615 

12055.015 

2090.090 

2909.950 

965.968 

912.368 

905.594 

905.534 

6 

26.892 

5.861 

4.092 

2.505 

.537 

.011 

.002 

0 


This method has been used in routine computation for a year 
without encountering a divergent case in a practical problem. At 
least for special cases when temperature is assigned, the process 
will converge for all values of the first estimates. Divergence 
is known to occur for certain cases where temperature is used as 
a variable when the first estimate of temperature and composition 
is sufficiently in error. Although no mathematical analysis has 
been made to determine the theoretical limits of convergence, the 
process appears to be satisfactory for practical' computation. 

Special treatment would be required if divergence is encoun- 
tered. Obtaining convergence should be possible by a sufficiently 
close new estimate of composition and temperature . This procedure 
is recommended when it is feasible but other procedures are pos- 
sible, depending on the individual case. 


Isentropic Expansion to Fixed Pressure 

The temperature and the composition of the products of reac- 
tion following an isentropic-expansion ratio of 20.4 at chemical 
equilibrium were also computed for the products of reaction of 
this example. The value of s Q is found from equation (28) by 
using the final values of each constituent of the adiabatic com- 
bustion and the absolute entropy values corresponding to the 
final combustion temperature. The calculated value of s 0 was 
763.476 calories per °K per mole. 
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First estimates . - The number of approximations necessary for 
a complete calculation can be considerably reduced If the initial 
estimate is based on previous experience. The final values of 
n^ and A determined for the combustion process of this e xamp le 
can therefore be the basis for this first estimate. 


Because the expansion ratio is 20.4, the four largest con- 
stituents can be estimated to be 1/20.4 of their combustion value. 


Pbf 3 = 



A = 


0.1304 

0.3503 

0.0867 

0.2544 

0.0815 


For convenience of presentation, the temperature was estimated 
to be 4000° ,K so that the values of table I could be used again. 
The remaining constituents can be estimated from the dissociation 
equations by setting log kj_ = 0. For example, p F would be 
determined with the assumed values of Pgj, and pg- from 
equation (6) and table I 


0 = log 0.3503 - log 0.0867 - log - 1-8944 
log Pj, = -0.45556 + 1.06198 - 1.8944 
= -1.28798 


Pp, = 0.0515 


Similarly, 
and Pj, 



can be estimated with the assumed values of 


log 0.1304 
log Pg = 


- log p - 3 log 0.0515 - 5.6953 
B 

-0.88472 + 3.86394 - 5.6953 
0.0019 
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If this procedure is followed for all the remaining constituents, 
the following list of first estimates can he made: 


p bf 3 

= 0.1304 

P 0H 

= 0.0150 

>b 2°3 

= 0.0078 

PtTT3 

= 0.3503 

-‘bf 


P BF 

= 0.0043 

% 

= 0.0269 


p =0 

BH 


■BO 




'H, 


H 2° 


0.0053 

p^ = 0.0867 
a • 

0 

p^ = 0 . 0019 

0.0029 

Pj, = 0.0515 

0.0009 

p =0.2544 
0 

0.0815 

T = 4000° E 


Correction equations . - When these estimates are substituted 
in equations (8) to (11) and (29), the parameters a, b, c, d, P, 
and s are calculated to be 1.9325, 5.6393, 9.7828, 4.3288, 0.9383, 
and 766.297, respectively. The correction equations (given in 
matrix form, in fig. 5) are then determined when these estimates 
and parameters are used in equations (19), (21), and (30). The 
remainder of this problem is continued in the same manner as the 
combustion calculation. The final solution is obtained in three 
trials . 


Isentropic Expansion to Mach Number of 1 

The temperature and the composition of the products of reac- 
tion following an isentropic expansion to the local velocity of 
sound by assuming chemical equilibrium was computed for the pro- 
ducts of reaction considered in this example. The value of s Q 
is the same as that found for the isentropic expansion to 1 
atmosphere. 

First estimate . - For simplicity, the same first estimates 
of 1 mole, 1, and 4000° E for n^. A, and T, respectively, 
were again made. 


O 
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Correction equations. - The matrix shown in figure 6 was con- 
structed from the application of the data of table I to equations 
(19), (21), (30), and (44). The first sixteen rows are identical 
to the corresponding rows of figure 4 because the same equations 
are used ; the next row is obtained from equation ' (30) . These 
seventeen equations are then solved to determine the values of 

d log A 
d log T 

(last row of fig. 6). As in the previous calculations, the 
resulting corrections are then applied to the first estimates of 
n^, A, and T by using equation (56) and the process is repeated 

until the assigned conditions are satisfied. 

Lewis Flight Propulsion Laboratory, 

national Advisory. Committee for Aeronautics, 

Cleveland, Ohio, September 7, 1949. 


) , which are used in equation (44) 
s 
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APPE&DU A 
SIMBOLS 

number of equivalent foimulas 

animation of each, atomic type over products of reaction 
per equivalent formula; with subscript, number of 
atoms of each element within chemical formula 

molar specific heat at constant pressure and standard 
conditions 

specific heat coefficient for matrix 

molar specific heat at constant volume and standard 
conditions 

molar internal energy at standard conditions 
Internal energy per equivalent formula 
molar free energy at .standard conditions 
molar enthalpy at 0° K and standard conditions 
molar enthalpy at standard conditions 
enthalpy per equivalent formula 
enthalpy coefficient for matrix 

sum of heat and kinetic energies per equivalent formula 
equilibrium constant 

ratio of equilibrium constant based on partial pressures 
to equilibrium constant based on free-energy change 
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M 


m 




3SL 


N 


P 

Q 


<L,r 

B 


s 

s' 

T 

t 

U 

u 

V 

v 


Mach, number 

molecular weight of equivalent formula 
mass flow per second 
number of dissociation equations 
number of chemical elements involved in reaction 
number of products of reaction 
‘.number of moles 
total pressure 
partial pressure 
any function 
any variables 
gas constant 

molar entropy at standard conditions 

entropy per equivalent formula 

entropy coefficient for matrix 

temperature 

throat area 

unit matrix 

velocity of sound 

volume 

velocity of flew 
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• • • f jLj X, Z 


x 

a 

5 

P 


elements within representative chemical formula 

matrix variable 

submatrix 

total -error parameter 
density 


Subscripts : 
a, b, c, • • • 
f 

g 

1 

m 

o 

P 

s 

T 

• • • y X, T y Z 


number of atoms within chemical formula 

fuel 

oxidant 

any point in nozzle 

number of types of gaseous molecule 

initial given condition 

constant pressure 

constant entropy 

temperature, °K 

product index numbers (i) that designate atomic 
gases 



i 


product index number 


f 
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APPENDIX B 


I MATRIX CONSTRUCTION AND REDUCTION 

I 

) 

A coefficient matrix is a scheme of detached coefficients of 
a set of linear equations that are to "be solved simultaneously. 

An augmented matrix is identical to a coefficient matrix except 
that the constants are included. Equations (19) to (23) con- 
stitute such a set of equations for the simultaneous determin- 
ation of the variables A log n^, A log A, and A log T. 

Construction. - Because of the large number of zeros occur- 
ring in the matrix, a considerable saving in effort can be made 
by proper arrangement of the order of the rows and the columns. 
The following arrangement provides a partly symmetrical matrix 
that has been found to be among the easiest to evaluate as long 
as the products of reaction are principally gaseous and the 
dissociation constants are expressed in terms of the atomic 
species : 

1. The order of the columns should be 

(a) A log n^ of gaseous molecules 

(b) A log n^ of atoms 

(c) A log of liquid and solid products 

(d) A log A 

(e) A log T 

(f) Constant terms of equations 

2. The order of the rows is 

(a) Dissociation equations in same order as gaseous 
molecules in columns 


(b) Mass-balance equations in order of atoms in columns 
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(c) Dissociation equations for solid and liquid pro- 

ducts in same order as solid and liquid in 
columns 

(d) Total-pressure equation 

(e) Heat-tall ance equation in combustion calculation; 

entropy-balance equation in calculation of 
isentropic expansion to fixed pressure 


In the calculation of isentropic expansion to an assigned 
Mach number, the order of the columns is not changed hut the 
rows are modified as follows: (l) The entropy-balance equation 

is substituted for the total -pressure equation; and (2) the 
heat-balance equation is changed to include the kinetic energy 
in accordance with equation (44). 

The -values of f - a i°8 A ) are first com- 

V d log T \ a log t/ b 

puted by means of the Hp+1 order matrix and are then substituted 
in equations (45) and (46) to yield the (Bp+2)nd row of the com- 
plete matrix. 


Solution. - One of the best methods of solving simultaneous 
linear equations is given by Crout (reference 8) . With this 
method, an auxiliary matrix is constructed from an original aug- 
mented matrix by a simple routine. T h i s auxiliary matrix is of 
the order equal to the original matrix. The solution for the set 
of equations can be obtained by a process of back substitution in 
the auxiliary matrix. 


For convenience, the order of the matrix is reduced before 
the Crout method is applied. A matrix arranged as recommended 
can be partitioned so that a unit matrix of the order (m,m) 

appears in the upper left corner, where m is equal to the number 
of types of gaseous molecule. The original augmented matrix can 
then be written 


i 

% ' a i 


a ' a 
2 i 3 


(Bl) 
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When the- Crout method is applied to the original augmented matrix, 
the Crout auxiliary matrix can he expressed as 


_ u ? i X _ 
a 2 : V 


(B2) 


■where uj , JaJ , and jjXgj are identical, to the corresponding 
Buhma trices of the original matrix. By observing the operations 
involved in the construction of the Crout auxiliary matrix, jjx^J 
is shown to he identical to the auxiliary matrix of the augmented 
matrix px^J defined by 

■ [°s] ■ [ a s] ■ [ a s] [°l] (K5) 

For computation, equation (B3) is written 



(B4) 


where Ujj- is a unit matrix of order equal to the number of col- 
umns of [a~] . The numerical solution is then obtained by carrying 


out the matrix multiplication indicated in equation (B4) to find 
Jogj . The Crout auxiliary matrix is constructed from Jagj 

The values of the variable x(m+l), . . x(H +2) are found 


from by the process of back substitution given by Crout. The 

values of the remaining variables are found by the matrix equation 


~ X 1 

• 


V 1 

• 

•* 

• 

jSa 

“'[“I] 

• 

• 

%+2 

_-l 
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For illustration, the submatrices Jjx^J , |~a^J , and 
were taken from figure 4 and used to construct figure 7 . The last 
six rows of figure 4 correspond to the matrix j^oc 2 jc 

shown in figure 7(a). The matrix 


- a l 


^m 


, ia^j and axe 


is shown in figure 7(b) 


where, for convenience of computation, the columns have been tab- 
ulated as rows with the first row at the top. The opera tians 

required by equation (B4) to evaluate . JbtjTJ are indicated in fig- 
ure 8(a) and the results for the example given in figure 7 are shown 

in figure 8(b). The operation ^ ( © ©■) is illustrated in detail 
as follows : 

0X0+0X0+0X0+1X 1.000 +0X0+0X0+2X 2.000 + 

2 X 2.000 + 1 X 1.000 + 1 X 1.000 + 0X0 + 0X0 + 

1 X 1.000 + 0x0 + 0x0 + 0x0 = 12.000 


Practical computation . - la practical computations, writing 
the complete original matrix, as shown in figure 4, 1 b unnecessary; 
instead, the matrices shown in figure 7(a) and 7(b) are written out 
so that figure 8(b) may be obtained. Except for the laBt two rows, 
figure 7(b) is always the same for a given group of chemical elements. 

The process of obtaining figure 7 is aB follows : 


1. Value 8 of are entered in row E. 

2 . Values of the elements of rows (A) , (B) , (c) , and (5) for 
columns 1 through N are obtained by multiplying the elements of 
row (E) by the elements of row (l), (?), (3), and (£), res- 
pectively. 

3. Values of the elements of row (f) for columns 1 through N 
are obtained by multiplying the elements of row (g) by the values 
of (ttji)jL found from a table of the thermodynamic properties for 
the substances. 
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4. Column A log A is found "by summing the elements in each 
row and writing the negative of the total in column A log A 
except for row ® where the value is 0. 

5. The elements of the A log T and constant col umn s are 
evaluated by means of the expressions shown in figure 1. 

6. Bow (D is obtained by entering the values of (AH^/RT from 
a table of thermodynamic properties. 

7. Elements of row (?) are obtained from the equations defining 
•log kj[, by taking log Kj. from a table of thermodynamic properties. 


1. Brinkley, Stuart B. , Jr.: Calculation of the Equilibrium Com- 

position of Systems of Many Constituents. Jour. Chem. Phys., 
vol. 15, no. 2, Eeb, 1947, pp. 107.-110. 

2. Erieger, F. J., and White, W. B.: A Simplified Method for Com- 

puting the Equilibrium Composition of Gaseous Systems. Jour. 
Chem. Phys., vol. 16, no. 4, April 1948, pp. 358-360. 

3. Huff, Vearl N. , and Calvert, Clyde S.: Charts for the Computa- 

tion of ■ Equilibrium Composition of Chemical Beactions in the 
Carbon-Hydrogen-Oxygen-Nitrogen System at Temperatures 
from 20000 to 5000° K. NACA TN 1653, 1948. 

4. McCa nn , W. J., rev. by L. B. Turner, and Emory A. Bauer: Ther- 

modynamic Charts for Intemal-Combustlon-Engine Fluids. 

NACA TN 1883, 1949. 

5. Scarborough, James B.: Numerical Mathematical Analysis. Johns 

Hopkins Press (Baltimore), 1930, pp. 187-190. 

. Both, W. A., und Borger, Erika: Zur Thermochemie des Bors. 

Berichte d. D. Chemischen Gesellschaft, Jahrg. 70, Nr. 2, 

Jan. 6, 1937. S. 48-54. 
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7. Anon.: Tables of Selected Yalues of Chemical Thermodynamic 

Properties. Hat. Bhr. Standards, Dec. 31, 1947. 

8. Crout, Prescott D. : A Short Method for Evaluating Determin- 

ants and Solving Systems of Linear Equations with Beal or 
Complex Yariables. AIEE Trans. (Suppl.), vol. 60, 1941, 
pp. 1235-1241. 
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TABLE I - VALUES OP CONSTANTS FOR REACTION OP DIBORANE WITH FLUORINE 

OXIDE (BgHg + 5PgO) 


Product 

Fixed 

Determined at estimated temperature 
of 4000° K 

i 

a i 

b i 

c i 

d i 

(Hg) ± 

(AHO/RT^ 

(s^i 

< c S>i 

log 

Equivalent 

formula 

0 

2 



6 

10 

5 






BP 3 

1 

1 

0 

3 

0 

72.172 

-62.0753 

105.951 

19.738 

5.6953 

b 2°3 

2 

2 

0 

0 

3 

233.435 

-80.5932 

116.760 

25.660 

5.1094 

BP 

3 

1 

0 

1 

0 

262.961 

-17.2884 

73.904 

8.905 

1.6342 

BH 

4 

1 

1 

0 

0 

356.994 

-8.3004 

61.412 

8.826 

-2.6110 

BO 

5 

1 

0 

0 

1 

252.739 

-18.1834 

69.620 

9.065 

1.0327 

B 2 

6 

2 

0 

0 

0 
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Figure 1. - General matrix for solution of oorreotion equations for adiabatic combustion and isentropio expansion to assigned 
pressure. Equations: (19), dissociation of gaseous woleoules; (21), mass balance} (20), dissociation' of solids or 
liquids; (22), pressure; (25), heat balance; (30), entropy balance. All blank spaces denote Eeros, 
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Figaro 2 . - General matrix for doterninatiari of derivatives used for calculating the local velocity of sound. 
Squatlonai (35), dissociation of gaseous molecules; (37), mass balance; (36), dissociation cf solid or 
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Figure 3. - General matrix for solution of correction equations for the process of isentropio 
expansion to assigned Mach number. Equations! (19), dissociation of gaseous molecules; 
(21), mass balance; (20), dissociation of solid or liquid moleoules; (30), entropy balance; 
(44), energy balance. All blank spaces denote seros. 
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Figure 7. - Breakdown of complete matrix of exacplo to facilitate calculation. An h i s-nk spaces denote zeros. 


& 


\ 

! 

f 

l 


f 

r 


F 


I 


f 


8921 


NACA TN 2113 







































































NACA TN 2113 


47 


A log A log ng 


011111111 

A log A 

A log T 




2(®@ ) 

2 ( ®® ) 

2 ( ©0 ) 

2(©® ) 

2(@® ) 

2((TXD ) 

2 ( dxD ) 

2 ( dXD ) 

2 ( dXD ) 

2((D® ) 

2(©® ) 

2 ( ®® ) 

2( CD© ) 

2(d)© ) 

2 ( (D© ) 

2 ( (D© ) 

2 ( eg© ) 

2 ( ©© ) 

[VJ 

© 

@ 

2 ( CD® ) 

2 ( dxg) ) 

2(®® ) 

2 ( ®© ) 

2 ( dX§) ) 

2( dXe) ) 

2 ( @® ) 

2( CD©) 

2(0© ) 

2 ( @© ) 

2 ( ®© ) 

2 ( ©© ) 

2 ( ©©) 

2(@®) 

2(@® ) 

2 ( dXD ) 

2(©@ ) 

2(0© ) 

2(®© ) 

2 ( ©® ) 

2 ( @© ) 

(a) Method of calculation of reduced matrix. 

12.000 

1.000 

1.000 

3.000 

-8.000 

“ -127.873 

1.391 

1.000 

13.000 

4.000 

7.000 

-9.000 

-283.010 

-16.322 

1.000 

4.000 

16.000 

0 

-8.000 

-240.597 

-13.564 

3.000 

7.000 

0 

17.000 

-9.000 

-333.400 

-17.383 

8.000 

9.000 

8.000 

9.000 

0 

-294.871 

-3.866 

885.360 

2873.620 

786.118 

1241.423 

-2734.615 

-45475.661 

-766.321 
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Figure 8. - Method of reduction of order of example matrix. 
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